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We report on an extensive investigation of the multiferroic compound TbMnOs using x-ray scat- 
tering techniques. Non-resonant x-ray magnetic scattering (NRXMS) was used to characterise the 
domain population of the single crystal used in our experiments. This revealed that the dominant 
domain is overwhelmingly A-type. The temperature dependence of the intensity and wavevector 
associated with the incommensurate magnetic order was found to be in good agreement with neu- 
tron scattering data. X-ray resonant scattering experiments were performed in the vicinity of the 
Mn K and Tb L3 edges in the high-temperature coUinear phase, the intermediate temperature cy- 
cloidal/ferroelectric phase, and the low-temperature phase. In the coUinear phase, where according 
to neutron diffraction only the Mn sublattice is ordered, resonant El — El satellites were found 
at the Mn K edge associated with A-type but also _F-type peaks. Detailed measurements of the 
azimuthal dependence of the F-type satellites (and their absence in the NRXMS experiments) leads 
us to conclude that they are most likely non-magnetic in origin. We suggest instead that they may 
be associated with an induced charge multipole. At the Tb L3 edge resonant A- and _F-type satel- 
lites were observed in the coUinear phase again associated with El — El events. These we attribute 
to a polarisation of the Tb 5d states by the ordering of the Mn sublattice. On cooling into the 
cycloidal/ferroelectric phase a new set of resonant satellites appear corresponding to C-type order. 
These appear at the Tb L3 edge only. In addition to a dominant El — El component in the a — tt' 
channel, a weaker component is found in the pre-edge with a ~ a' polarization and displaced by 
—7 eV with respect to the El — El component. Comprehensive calculations of the x-ray scattering 
cross-section were performed using the FDMNES code. These calculations show that the unro- 
tated a — a' component of the Tb I/3 C-type peaks appearing in the ferroelectric phase contains a 
contribution from a multipole that is odd with respect to both space and time, known in various 
contexts as the anapole. Our experiments thus provide tentative evidence for the existence of a novel 
type of anapolar order parameter in the rare-earth manganite class of mulitferroic compounds. 



A. Introduction 

Materials possessing more than one "ferroic" order 
(which includes magnetism, ferroelectricity, elastic dis- 
tortions, etc.) are referred to as multiferrroic. Although 
multiferroics are rare — magnetism requiring a partially 
filled electronic shell, while canonical ferroelectrics such 
as BaTiOa have an empty one — they are of considerable 
interest both from a fundamental point of view, and from 
the possibilities that they offer in the field of spintronics. 
Recently the RMnOs and RMn2 05 (R=rare earth) se- 
ries of compounds have attracted a considerable amount 
of attention, as it has been shown that members of these 
series are multiferroic and display novel effectsii^i^i^i^i^. 
In the case of TbMnOs, for example, a giant magne- 
toelectric effect has been reported allowing the electric 
polarization to be switched by applied magnetic fields^. 
Correspondingly, it has been shown in HoMnOa that the 
magnetic response can be controlled by an applied elec- 
tric field^. Of key importance in understanding the func- 
tionality of these materials is to obtain a complete de- 
scription of the various order parameters. 

To shed further light on the nature of the order pa- 
rameters in TbMnOa we have employed x-ray scattering 



techniques, with an emphasis on exploiting the rich possi- 
bilities offered by x-ray resonant scattering (XRS). In an 
XRS experiment, the photon energy is tuned close to an 
absorption edge and an excited intermediate state is cre- 
ated by transitions between core and valence shells. At 
resonance, the scattering amplitude of a single ion can 
be described by a tensor, whose properties are related 
directly to specific terms in the multipolar expansions of 
the electric charge and magnetization^. XRS is now well 
established as being capable of revealing a great diver- 
sity of ordering phenomena in solids. These are usually 
probed via the pure excitation channels El — El and 
E2 — E2, where the two distinct channels arise from se- 
lection rules for the change in angular momentum AL 
between the core and excited states AL{E1 — El) — ±1 
and AL{E2 — E2) — ±2, respectively. For such events, 
the tensors are even parity and are ranked according 
to their behaviour under time reversal as time-odd, for 
magnetic, and time-even for charge orderings. If, how- 
ever, the ions are located in crystallographic positions 
which lack a centre of inversion symmetry (odd parity), 
hybridisation will occur between the valence orbitals of 
that atom and transitions may then occur via mixed pro- 
cesses to the hybridised states, giving rise to parity odd 
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FIG. 1; (Color online) a The crystal structure of TbMnOa in 
its high-temperature (Pbnm) phase shown in projection down 
the c axis. 



El — E2 events. Such events open up new and interest- 
ing possibilities for studying ordering phenomena driven 
by parity odd multipoles. For example, magnetoelectric 
toroidal moments may develop which can be visualised 
as arising from placing a series of magnetic dipole mo- 
ments end to end so as to form a closed circular pathS. 
Such moments are characterised as having odd symmetry 
with respect to parity and can have odd or even symme- 
try in time reversal. These toroidal moments have his- 
torically been referred to as anapoles first considered in 
the context of multipolar expansions in nuclear physics^. 
Multipoles associated with El — E2 processes have been 
shown to be of importance in non-reciprocal and natural 
circular dichroisroiSiiiiii^, and recent predictions having 
been made for XRS experiment a^'^'^'^i^^i^^ . 

At room temperature TbMnOa adopts an orthorhom- 
bic structure, space group Pbnm, with lattice parameters 
a=5.316 A, b=5.831 A and c=7.375 A (Fig.[J)^ The Mn 
and Tb ions occupy the 46 (point symmetry 1, i.e. in- 
version symmetry) and 4c (point symmetry m) sites, re- 
spectively. Neutron diffraction experiments have shown 
that the Mn'^+ magnetic moments order below T7vi=41 
K to form an incommensurate, longitudinally polarized 
structure (hereafter referred to as the coUinear phase) 
described by a wave-vector (0,qMn,0) with qAf„=0.28- 
0.29 b*— li^. In this coUinear phase the Tb moments are 
thought to be disorderedii. Below Tn2 = 28 K the mag- 
netic structure becomes non-collinear (hereafter referred 
to as the cycloidal phase) by developing a component 
of the magnetization on the Mn sublattice along the c 
axis, and at exactly the same temperature the material 
becomes ferroelectrioii. According to modelling of the 
neutron diffraction, the Tb moments in this phase or- 
der at the same wavevector as the Mn moments, but are 
transversely polarized along the a axis. On further cool- 
ing below Tn3 = 7 K, the Tb moments undergo a further 
transition and develop a component of the magnetization 
at a distinct wave- vector of (0, qxb, 0) with qrb=0.42 b*, 
again thought to be transversely polarized along the a 
axis. Kenzelmann et alr^ have argued that the transition 
to a ferroelectric state results from a loss of a centre of in- 
version symmetry when the magnetic structure changes 



from being coUinear to non-collinear at Tn2 = 28 K. 
This idea has generated considerable interest, and fer- 
roelectric transitions driven by the formation of non- 
collinear magnetic structures have now been reported in 
a number of systemsi^i^S. Several different theoretical 
approaches have been developed to explain the magneto- 
electric properties of this class of multiferroio2ii2^, which 
suggest that the macroscopic electric polarization P is 
related to the magnetic moment Si by 

P = aetj X {Si X Sj) 

where is a vector connecting the spins on the sites i 
and j, and a is a constant. 

It should be noted that a certain degree of inconsis- 
tency exists in the neutron diffraction literature con- 
cerning the magnetic structure adopted by TbMnOs. 
In TbMnOs the magnetic reflections may be classified 
into four distinct types, which have been shown to 
have distinct temperature dependencesiiiiSii^i^i. The F- 
type structure gives rise to magnetic reflections classi- 
fied hy h + /c=cven, ^=:even; the G-type reflections with 
h + A:=odd, Z=odd; the C-type h + A:=odd and ^=even; 
and the ^-type structure gives reflections corresponding 
to h + k—even, Z=odd. In the original investigation by 
Quezel et al^ strong ^-type and weak G-type magnetic 
peaks were reported, with the systematic absence of G- 
type and F-type reflections. Blasco et al.~ obtained 
similar results from powder neutron diffraction. In their 
neutron work on a large single crystal of TbMnOa, Kaji- 
moto et al}^ found evidence for all F-, C- and G-type 
structures. The A and G-type structures have their mo- 
ments unambiguously along the b axis, but attempts to 
determine the moment directions for the -F and G-type 
structures could only constrain the direction to be either 
the c or a axis. The more recent neutron investigation by 
Kenzelmann et al.^ reported only the ^-type reflections 
with moments along the b axis in the coUinear phase, in 
agreement with Kajimoto et al. A plausible explanation 
for the observation of the different magnetic structures in 
these experiments is that they depend on details of sam- 
ple growth and preparation. All neutron investigations, 
however, agree on one thing: the ^-type structure is al- 
ways the dominant phase. It is also worth noting that 
the magnetic order is accompanied by weak charge satel- 
lites at 2qMn which have been observed with neutron and 
non-resonant x-ray scattering22i2^. 

Our experiments on TbMnOa were undertaken in the 
spirit of exploring whether the well documented charac- 
teristics of XRS could be exploited to reveal information 
complementary to that provided by neutron diffraction. 
In particular, the element and electron shell specificity of 
x-ray resonant magnetic scattering (XRMS) is often used 
to isolate the contribution from individual components 
in systems with more than one type of magnetic species, 
such as TbMnOs. A second major objective of our study 
was to investigate whether the loss of inversion symmetry 
that must accompany the transition to the ferroelectric 
state in TbMnOs allows new terms in the XRS cross- 
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FIG. 2: (Color online) Definitions of the nomenclature used 
to label the polarization state of the photon beam and the 
coordinate system used to resolve the components of the mag- 
netic moments relative to the incident and scattered beams. 



section to become visible, such as El — E2 terms associ- 
ated with the ordering of parity odd multipole a^^'^"^i^^i^^ . 

This paper is organised as follows. In the next sec- 
tion, [21 we outhne the experimental details, including 
a description of the scattering geometries (crystal ori- 
entation, x-ray polarisation, definition of azimuth, etc.) 
employed in our experiments. In Sec. [C] we present the 
experimental results, which for clarity of exposition is 
further sub-divided into sub-sections. Our attempts to 
characterise the magnetic domain structure of our sample 
using non-resonant x-ray magnetic scattering (NRXMS) 
are described in Sec. I C II The results of the XRS in- 
vestigations in the coUinear phase at the Mn K-edge are 
described in Sec. I C 2i while in Sec. I C 31 we describe the 
outcome of XRS experiments in the coUinear phase at 
the Tb L3 edge. XRS investigations undertaken in the 
cycloid phase, again at the Mn K and Tb L3 edges, are 
described in Sec. I C 41 and Sec. I C 5[ respectively. Results 
taken in the the low temperature regime, T < 8K, at 
the Tb L3 edge are presented in Sec. I C 61 The inves- 
tigation of the 2q charge satellites is presented in Sec. 
I C 71 Our results are analysed and discussed in Sec. [Dj 
in the light of calculations of the electronic structure us- 
ing the FDMNES code^^. In the final section, [El we 
summarise our work. 



B. Experimental Details 

Single crystals of TbMnOa were prepared in Oxford by 
the floating zone method in an image furnace. The crys- 



tals were cut and polished to produce samples with either 
the b or a axis as the surface normal. The samples were 
investigated by SQUID magnetometery and were found 
to display the same behaviour as previously reported'^S. 
The lattice parameters refined in our synchrotron exper- 
iments were in very good agreement with with those de- 
duced from neutron scattering results and the our crystal 
mosaic was found to be w 0.03 degrees. These observa- 
tions indicate that the samples used in our studies were of 
execellent quality with comparable bulk magnetic prop- 
erties to other samples described in literature. 

The scattering experiments were undertaken at the 
XMaS and ID20 magnetic scattering beamlines at the 
ESRF, Grenoble, France. XMaS views radiation from a 
bending magnet and employs a vertical scattering geom- 
etry (cr-polarised incident photons). The sample used on 
XMaS was cut and polished with a b axis as surface nor- 
mal (see Fig. [5] for the definition of photon polarizations 
and coordinate system used in our experiment). The 
experiments on the undulator beamline ID20 were per- 
formed also mainly with the same vertical geometry as 
on XMaS. Some additional experiments on ID20 were 
carried out in a normal beam geometry, with the scatter- 
ing plane approximately horizontal, corresponding to tt- 
polarised incident photons. XRS measurements were un- 
dertaken around the Mn K (6.552 keV) and Tb L3 (7.515 
keV) edges, using a Cu(220) or Au(222) polarisation 
analyser, respectively. The wave-vector resolution was of 
the same order of magnitude in all our x-ray experiments, 
typically « 1x10^"* reciprocal lattice units (r.l.u.). The 
photon flux on ID20 is approximately 1x10^^ per sec- 
ond in a focused spot size of about 0.3x0.4mm^ (verti- 
cal X horizontal); on XMaS the corresponding numbers 
are 5x10^"'^ per second in 0.8x0.5mm^ at x-ray energy 
of 8keV and 200 mA current in the ESRF storage ring. 
The intense photon flux on ID20 carries with it the po- 
tential to heat the sample. Checks were made for such 
an effect, and when present (typically for sample tem- 
peratures below 10 K), the incident beam was attenu- 
ated by up to an order of magnitude. A key aspect of 
our experiments has been the utilisation of the high flux 
from the ID20 undulator source, which enabled us to per- 
form non-resonant x-ray magnetic scattering (NRXMS) 
to characterise the magnetic domains present in the near 
surface volume probed by the x-rays. These experiments 
were undertaken with incident energy of 7.470 keV using 
a Au(222)polarisation analyser. 

In the following, we use the conventional notation to 
describe the polarisation conditions of the incident and 
scattered photon beams, a-a' and tt-tt' for incident po- 
larisation scattered without rotation and ct-tt' and tt-ct' 
for photons scattered with a 90°-rotation of their electric 
vector. 

In our experiments extensive use was made of scans of 
the azimuthal angle "0, i.e. rotation around the scatter- 
ing vector—. In recent years the azimuthal dependence of 
XRS has been exploited increasingly in studies of systems 
displaying magnetic and multipolar order. The reason 
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FIG. 3: (Color online) Summary of azimuthal scans of vari- 
ous diffraction peaks recorded at different temperatures and 
photon energies. The data were recorded in the a — n' channel 
and have been normalised to unity. The lines are guides to 
the eye. 



is that the variation of the resonant intensity in a scan 
of the azimuthal angle reflects directly the symmetry of 
the intermediate states (orbitals) through which the reso- 
nance proceeds. By way of example, we show in Fig.[3I^a) 
the azimuthal dependence of the space-group forbidden 
(3,0,0) reflection recorded at the Tb L3 edge on XMaS. 
The intensity displays a cos^ dependence in agreement 
with calculations of E\ — El Templeton scattering^^!^ 
for the Pbnm space group. In either the experiments on 
samples with b or a faces as the surface normal direction, 
the crystal was oriented with the c axis in the scatter- 
ing plane, which defined the ^=0 of the azimuthal scans. 
Due to the construction of the diffractometers used in our 
experiments, azimuthal scans could only be performed in 
the vertical scattering geometry, and then mainly for re- 
flections parallel to the b axis. 



C. Experimental Results 

1. Magnetic Structure Investigation with NRXMS 

In order to ascertain which magnetic structures were 
present in our sample, and in particular to under- 
stand any effects that may arise in the near surface re- 
gion probed in our x-ray measurements, we initially in- 
vestigated the non-resonant x-ray magnetic scattering 
(NRXMS). Time constraints, and the weakness of the 
signal, prevented us from making exhaustive studies in 
the NRXMS regime. Nevertheless, these investigations 
turned out to be helpful in the interpretation of our XRS 
results. 

The NRXMS experiments were undertaken at a tem- 
perature of T—IA K where all the possible A-, F-, C- and 
G-type reflections were found to be present in the study 
of Kajimoto et al.^"^. The crystal was orientated initially 
at an azimuth of ip = 0, and a photon energy of 7.470 
keV was selected, below the Tb L3 edge. The scattered 
intensity calculated for NRXMS can be written as^ 



(1) 



where Li, Si, S2 and S3 are the components of orbital 
and spin moments along orthogonal axes Ui, U2 and U3 
(Fig. [2|), is the scattering Bragg angle, and Yd repre- 
sents the domain volumes of the d = A, F, C, 01 G- 
type structures within the illuminated sample volume. 
Upon azimuthal rotation, consideration should be given 
in Eq. [T] for the change in the magnetisation compo- 
nents along the Ui, U2 and U3 axes, so that for example, 
S'2(i/)) = S2{^ = 0) cos -Sli^ = 0) sini/.. 

The results of our NRXMS investigations are sum- 
marised in Fig. m and Fig. [5] for the a-n' and a-a' po- 
larised intensities, respectively. From Eq. [Tl it can be 
seen that the (t-tt' channel is sensitive to magnetic mo- 
ments in the scattering plane along Ui and U3, and the a- 
a' channel is sensitive to magnetic moments directed per- 
pendicular to the scattering plane along C/2. In Fig. 21 sig- 
nificant NRXMS intensity is evident for wavevectors cor- 
responding to the A-type magnetic structure, (0, 4±q, 1) 
and (0, 6 ± q, 1). (The slightly weaker intensities found 
at (0, 6 ± q, 1) compared to (0, 4 ± q, 1) is 

probably due to the decreasing magnetic form-factor 
rather than the geometric terms in Eq. [1]) Scans, taken 
at wavevectors corresponding to the -F-type and C-type 
structures, at (0, 4 ± q, 0) and (0, 3 ± q, 0) are shown in 
Fig. Hl^c) and (d) , where no peaks are visible above back- 
ground. Turning to the a-a' channel, summarised in Fig. 
[5l we found no evidence for NRXMS for the wavevectors 
investigated, indicating that the magnetic moments in 
this phase are confined to the b — c plane. 

Thus the fact that we observe satellite reflections in 
our NRXMS experiments at positions (h,k ± q, /), with 
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FIG. 4: (Color online) Summary of the non-resonant x-ray 
magnetic scattering observed on ID20 in scans of wavevector 
transfer parallel to b* through various positions in reciprocal 
space associated with the A-, F-, or C-type magnetic struc- 
tures. The data were recorded in the ferroelectric/cycloidal 
phase at r = 14K, at an azimuth of = and in the a — tv' 
channel. 



h + k even, / even, appears to indicate that our sample is 
predominantly A-type. While the fact that we do not ob- 
serve any peaks in the unrotated a — a' channel leads us 
to conclude that in our sample, and at a temperature well 
inside the cycloidal phase, the magnetic moments are in 
the b — c plane. In this way our sample appears to share a 
similar domain configuration to the one studied by Ken- 
zelmann et al. who found that in this phase the moments 
form a cycloid in the b — c plane. Further information 
on the orientation of magnetic moments in our sample 
was provided by studying the azimuthal dependence of 
the (0, 4 ± q, 1) satellite, shown in Fig. [3)Jc). The inten- 
sity depends only weakly on azimuth, which from Eq. 
[1] implies that the moments have their strongest com- 
ponent along the b direction. It should be noted that 
the weak azimuth dependence observed could also arise 
from 2 magnetic domains oriented along both the a and 
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FIG. 5: (Color online) Summary of the non-resonant x-ray 
magnetic scattering observed on ID20 in scans of wavevector 
transfer parallel to b* through various positions in reciprocal 
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c axes in our scattering geometry. However, the absence 
of magnetic scattering in cr — ct' channel is not consistent 
with such a scenario. 

The temperature dependence of the intensity of 
NRXMS observed at (0,4 ± q, 1) is shown in Fig. [^a), 
where it appears to agree with that deduced from neu- 
tron scattering experiments, increasing steadily below 
7^1=41 K. We also find good agreement with neutron 
scattering for the thermal evolution of the magnetic 
wave- vector as shown in Fig. El^b), starting at a value 
of about 0.29 at T^vi, decreasing in value as the sample 
is cooled through the coUinear phase, before appearing 
to lock into a value of about 0.28 below Tn2 in the ferro- 
electric phase. However, in agreement with the neutron 
data, it is clear that this is not a genuine lock-in transi- 
tion, as below Tt^2 the value of the wavevector increases 
slightly down to base temperature. 

The conclusion to be drawn from our NRXMS inves- 
tigation is that the dominant magnetic structure in our 
sample appears to be A-type, with the magnetic moment 
oriented mainly along the b axis. No evidence is found 
for the F-type or C-type structures, presumably because 
either they are not present, or their domain volumes are 
too small to allow them to be observed. 
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all recorded in a — tt'. For ease of comparison, the intensities 
have been normalised to unity at low temperature. The data 
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2. Collinear phase: Mn K-edge 

Turning to our XRS studies, we first investigated the 
response with incident energy tuned close to the maxi- 
mum of the Mn fluorescence spectrum (Fig. [7{a)), and 



with the sample cooled into the collinear phase (Tj 



N2 



< 



T < Tjvi). Under these conditions, scans in reciprocal 
space were performed to search for the existence of satel- 
lite peaks. Well defined diffraction satellites were found 
(Fig. [SJc) and (d)) at nominal A- type magnetic wave- 
vectors such as (0, fc ± q, 1), with q « 0.285 b*, close to 
the value of qMn deduced from neutron diffraction exper- 
iments. An additional set of satellites were discovered at 
F-type positions, (0, A:±q, 0) k even (Fig. [5]^a) and (b)). 
Both sets of satellites were found in the rotated a — tt' 
channel only. Extensive searches at C-type positions such 
as (0, 3 ± q, 0), failed to find any evidence for XRS at the 
Mn K edge in this phase. The observation of both A- and 
F-type reflections at the Mn K edge is difficult to recon- 
cile with the NRXMS presented in Sec. I C ll where peaks 
attributable to the A- type structure only were present. 
Moreover, the data in Fig. [5] also show similar intensity 
at the A- and F-type reflections in the collinear phase. 
This would appear to undermine any claim that the lat- 
ter is magnetic scattering from a minority F-type phase: 
neutron scattering studies always show that A-type is 
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ity of the Mn K edge, (b)-(c) Representative energy depen- 
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ometries. The data were taken at ID 20. 



dominant. We comment on possible explanations for the 
discrepancy between the NRXMS and XRS later. 

In Fig. [T^b) we show the energy dependence of the 
(0, 4 - q, 0) at 30 K, for incident a and the (4, -q, 0) for 
TT-polarised photons (Fig. [Tl^c)), as the photon energy is 
tuned through the Mn K edge. The resonant response 
is seen to exhibit a maximum at an energy of 6.557 keV, 
slightly above the position of the flrst derivative in the 
absorption at 6.552 kcV. The energy of the resonance 
allows us to identify unambiguously the maximum with 
El— El dipole transitions, which at the Mn K edge probe 
the 4p states. (A weak, pre-edge feature at 6.540 keV is 
also evident.) The scattering was found to be mostly 
in the rotated channel. The scattering in the unrotated 
channel is more difficult to interpret with any certainty, 
due to changes in the diffuse charge scattering and the 
onset of fluorescence at the edge. 

To further address the nature of the order revealed by 
the F-type peaks, the azimuthal dependence of the scat- 
tering of the (0, 4— q, 0) at the Mn K edge was determined 
with the results shown in Fig. ^h) for the ct-tt' channel. 
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FIG. 8: (Color online) Representative scans of the wavevector 
transfer with the photon energy tuned to the peak of the 
resonance at the Mn K edge. The scans shown are through 
the positions of F-type ((a) and (b)) and yl-type ((c) and (d)) 
reflections, and were recorded in the a — tt' channel at ID20. 



FIG. 9: (Color online) Representative scans of the wavevector 
transfer with the photon energy tuned to the peak of the 
resonance at the Tb L3 edge. The scans shown are through 
the positions of f-type ((a) and (b)) and A-type ((c) and (d)) 
reflections, and were recorded in the a — tt' channel at 1D20. 



Tlie scattering is seen to display tlie same cos^ ip depen- 
dence as the (3,0,0) space group forbidden reflection (Fig. 
ini;a)). The El - El XRMS amplitude can be expressed 
as^ 

j^Mag / cr — (t' tt — cr' \ 

res \^ (T — tt' tt — tt' J 

/ zi cos ^-t-Z3 sin ^ \ 

\ Z3 sin f?— Zi cos ^ — Z2 sin 2^? / \ / 

Here the amplitude is written as a matrix in a basis of 
the linear components of the polarization perpendicular 
(ct) and parallel (tt) to the scattering plane, for the in- 
cident (unprimed) and scattered (primed) beams and z 
is a unit vector parallel to the magnetic moment, with 
Cartesian components defined with respect to the or- 
thogonal axes J7i, U2 and U3, (Fig. [2]). The above equa- 
tion establishes that El — El XRMS is forbidden in the 
cr — cr' channel, and is sensitive to the magnetic moments 
in the scattering plane for the ct-tt' channel. Thus the 
two-fold azimuthal symmetry we observe at (0, 4 — q, 0) 
(Fig. (Sib)), with the maximum at -0 = 0, could be con- 
sistent with magnetic scattering from an i^-type struc- 
ture for which the magnetic moments would have to 
be orientated along the c-axis. This moment direction 
is in agreement with that deduced from neutron scat- 
tering experimentsiS. However, when we compare our 
expectations for the NRXMS response that would arise 
from such an i^-type structure, with moments orientated 
along the c axis, we fail to find consistency with our own 
NRXMS results for the absence of signal at (0, 4 -I- q, 0). 
From Eq. [21 using the Bragg angles at the Mn K-edge 



0(0,4 - q,0)=40 and (0,4 - q, 1)=29.7, the measured 
XRMS intensity ratio I((0,4 q, 1))/I((0,4 q,0)) can 
be used to obtain the expected volume ratio Va/Vf for 
the F-type and A-type structures. This leads to the ra- 
tio V^/Vf=3. We can then plug this ratio into Eq. [1] 
for the non-resonant scattering described in Sec. I C 1[ 
using the NRXMS Bragg angles (0,4 - q, 0)=34.8 and 
(0, 4-q, 1)=23.8 and calculate the NRXMS intensity we 
would expect at (0, 4 4- q, 0). However, from this we find 
that the NRXMS should be significantly larger than that 
observed at (0, 4 -|- q, 1), not zero as we have found. 

Further doubt that the XRS at (0, 4±q, 0) is magnetic 
in origin comes from its temperature dependence. The 
thermal evolution of the (0, 4 — q, 0) ct-tt' intensities are 
shown in Fig.[6Ua) and appear to follow closely the A-type 
(0, 4 — q, 1) magnetic intensities that we measured using 
NRXMS. The close agreement in the temperature depen- 
dence is in contrast to that revealed in neutron diffraction 
experiments^^. In these measurements the temperature 
dependence for the two structures is very different: the 
F-type intensities increase at a slow rate below Tjvi=41 
K and then increase dramatically below Tjva = 7 K due 
to the apparent stabilisation of the F-, G- and C-type 
structures by the onset of Tb magnetic order. 

Also compared in Fig. [Hfb) is the temperature depen- 
dence of the modulation wavevector at (0, 4 — q, 0) with 
the non-resonant scattering at (0,4 — q, 1). Of signifi- 
cance, especially in the light of evidence that the XRS 
observed at (0, 4 — q, 0) may not be magnetic in origin, 
is the similar thermal evolution of wave-vector to that 
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FIG. 10: (Color online) (a) Fluorescence spectrum in the 
vicinity of the Tb L3 edges, (b)-(c) Representative energy 
dependence of the x-ray resonant scattering near the Tb L3 
edges at various satellite positions, and for various polariza- 
tion geometries. The data were taken on XMaS. 



found for A-type magnetic scattering we measured at 
(0,4 - q, 1) with NRXMS. This effect may suggest that 
the underlying order parameter is driven by the onset 
of magnetism of the Mn ions, suggesting an interesting 
coupHng of magnetic and charge order parameters in this 
material. 

In summary, in the coUinear phase at the Mn K edge 
we find ^-type reflections at (0, 4 — q, 1), which was the 
dominant peak in the NRXMS study. Its dependence 
on energy, polarization, temperature, etc., have the hall- 
marks of an XRMS process, where the ordering of the 
local Mn 3d moments leads to a polarization of the Ap 
band. The peak at the F-type position (0,4 ± q, 0) is 
less easy to interpret. If it is attributed to also being 
an XRMS peak, then we can infer that for the i^-type 
structure the Mn moments are polarized along the c axis. 
However, the characteristics of the (0, 4±q, 0) — not least 
of which is its absence in our NRXMS study — forces us 
to consider other possible origins for this peak. If this 
order parameter is not magnetic order, and given that it 



can be associated with El — El events, we may suppose 
that the XRS arises from El — El time-even events as- 
sociated with a charge multipole. Moreover, given that 
its wavevector follows so closely with the magnetic, this 
order parameter would have to be coupled linearly to the 
magnetism of the Mn sublattice. 



3. CoUinear phase: Tb L3 edge 

Experiments were also performed in the coUinear phase 
for photon energies in the vicinity of the Tb L3 edge. 
Resonant satellites were found at positions correspond- 
ing to both the nominal Mn A-type magnetic satellites 
at (0, 4 ± q, 1) and the seemingly anomalous peaks at 
(0, 4±q, 0). Representative scans of the wavevector trans- 
fer parallel to b* are shown in Fig. [51 where sharp diffrac- 
tion peaks are evident reflecting long-range nature of the 
ordering probed. Polarization analysis of these showed 
that the scattering was predominately in the rotated ct-tt' 
channels. The observation of these satellites is particu- 
larly notable as, according to all previous neutron scat- 
tering studies, the Tb moments in the coUinear phase are 
disordered. The data in Fig. [DJa) and (b) also appear to 
show that far from being the weaker peaks, at the Tb L3 
edge the i^-type are actually considerably more intense 
than the A-type peaks. Such comparisons should, how- 
ever, be qualified by the fact that at the Tb L3 edge a 
marked asymmetry in the intensity of the ±q satellites 
was present (see Fig. [TTTb')'). 

The energy dependence of (0, 4 ± q, 0) is shown in Fig. 
[rUf bV Inspection of the energy scans allows us to ascribe 
the resonances to an i?l — i?l processes, which in the case 
of the Tb Ls edge connect 2p to 5d states. The azimuthal 
dependence of the (0, 4 ± q, 0) ct-tt' intensities is shown 
in Fig. [TlTb) , and can be seen to exhibit the same cos^ t/j 
dependence as found at the Mn K edge (Fig.[3jb)). The 
temperature dependence of the (0, 4±q, 0) intensities and 
wave- vectors are shown in Fig. [T^] where we find that they 
follow exactly the temperature dependence of the XRS 
observed at the Mn K edge. 

From these observations, we can immediately deduce 
that in the interval T/v2 < ^ < T/vi, in other words well 
above any ordering of the Tb'^^ dipole moments, the 5d 
states at the Tb sites are polarized by ordering of the Mn 
magnetic moments and with a modulation wave-vector 
identical to that of the Mn sublattice. (In this context 
we note that induced polarisations have been reported on 
the anion species of several actinide materials"^^.) This of 
course does not answer the question of the nature of the 
polarization that is created at the Tb sites, and hence 
the origin of the scattering process which gives rise to 
the resonances. We consider the various possibilities in 
Sec.[Dl 
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FIG. 11: (Color online) Summary of azimuthal scans of var- 
ious diffraction peaks recorded at different temperatures and 
photon energies. The data were recorded in the a — n' channel 
and have been normalised to unity. The lines are guides to 
the eye. The azimuth t/j = corresponds to the c axis in the 
scattering plane. 



4. Ferroelectric/ cyclotdal phase: Mn K 



It is of considerable interest to follow the XRS at the 
Mn K edge on entering the ferroelectric phase. In the 
model proposed by Kenzelmann et al., ferroelectricity 
arises when the Mn sublattice makes the transition from 
a collinear to non-collinear structure, removing the in- 
version centre at the Mn site. In principle, the loss of 
inversion symmetry may be detected directly in XRS ex- 
periments as the opening of new E1 — E2 scattering chan- 
nels in the region of the pre-edge^^^ii^ii^. 

The most complete data set that we have on changes 
to the XRS spectra at the Mn K edge on entering the 
ferroelectric phase was taken with incident tt polarized 
photons. The temperature dependence of the resonant 
line shape of the i^-type satellite (4, — q, 0) is summarised 
in Fig. [131 In the collinear phase at 30 K the scattering 
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FIG. 12: (Color online) Temperature dependence of the 
intensity and modulation wave-vector at the Tb I/3 edge 
recorded in a — tt' . For ease of comparison, the intensities 
have been normalised to unity at low temperature. 



is predominantly in the rotated n — a' channel, with a 
smooth, weaker response in the unrotated channel, most 
likely associated with the onset of fluorescence at the 
edge. On cooling into the ferroelectric phase a dramatic 
change occurs in the unrotated channel with a peak ap- 
pearing displaced to lower energy by 2 eV from the peak 
in the response in the rotated channel. With further cool- 
ing the intensity of the unrotated component increases 
until it becomes comparable in magnitude to the rotated 
one. In addition a feature appears above the edge in both 
channels. 

Inspection of Eq. [5] reveals that for n polarized inci- 
dent photons. El — El XRMS scattering from a longi- 
tudinally modulated magnetic structure polarized along 
b should occur in the tt — tt' channel only. This is ex- 
actly the opposite of what is shown in Fig. [T51a). where 
the predominant scattering occurs in the tt — a' chan- 
nel. Further consideration of Eq. [2] indicates that, if the 
peak is to be ascribed to XRMS, then the moments in the 
collinear phase are instead polarized along a or c. Consis- 
tency with the data on the (0,4±q,0) satellites (Sec. I C 2"|) 
can be achieved when we recall that azimuthal scans at 
this peak position suggested moments polarized along c. 
The appearance of significant intensity in the unrotated 
channel on entering the ferroelectric phase (Fig. [TST b)) 
might then further be explained, for an El — El XRMS 
process, by the development of a component of the mag- 
netization along the b axis. However, this interpretation 
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FIG. 14: (Color online) Energy dependence of XRS at the 
(0,3-qAfn,0) satellite in the vicinity of the Tb L3 edge recorded 
on ID20. The data were taken at a temperature of r=14 K, 
and the azimuth was set to he — 20°, i.e. close to the 
minimum in the intensity of the response in the cr — vr' channel 

(Fig.mrb)). 



40 



20 



(c) (4 -q^_^,0) Mn K 
T=4 K I*- 



i 5 i 



7t-0 

Tt-n' 



6.52 



6.53 



6.54 6.55 6.56 
Energy [keV] 



6.57 



6.58 



FIG. 13: (Color online) Energy dependence of XRS at the 
(4,-qMn,0) satellite in the vicinity of the Mn K edge collected 
using a normal beam geometry with tt incident photons on 
ID20. 



would seem to be untenable as the intensity in the unro- 
tated channel peaks 2 eV below that in the rotated one, 
indicating that the former is not pure El — El. Here 
the caveats raised in Sec. I C 21 and I C 31 about interpret- 
ing the J^-type peaks as XRMS should also be recalled. 
We do not observe any C-type XRS peaks at positions 
such as (0,3±q,0) at the Mn K edge, which indicates that 
these Mn magnetic domains are absent in our sample, in 
agreement with our NRXMS results. 



5. Ferroelectric/ cycloidal phase: Tb Lz edge 

Upon cooling into the cycloidal phase Tn2 <T< Tni, 
a new set of C-type resonant satellites at (0, fc ± q, 0) k 
odd were revealed at the Tb edge. These were found 
to be absent at the same temperature at the Mn K edge 
and unlike the A-type and i^-type reflections, we do not 
find any charge order satellites at 2qMn associated with 
these C-type peaks. Key data on their energy depen- 
dence measured at (0, 3-l-q, 0) are shown in Fig. [HI where 
it is apparent that the intensity appears mostly in the 
rotated a-ir' channel. The azimuthal dependence of the 
(0,3 + q, 0) is shown in Fig. [TlTb) and Fig. [THl where 
it is seen to have a sin^ ip dependence, in anti-phase to 
the azimuthal dependence found at (0, 4 -f q, 0) shown in 
Fig. Elb) and Fig. [TTTa). At these C-type positions, en- 
ergy scans also reveal a weaker pre-edge feature as seen 
in Fig.[T31 which is scattered with a-a' polarisation. The 
azimuthal dependence of this peak has a cos^ ip azimuthal 
dependence (Fig. [T5)) . 

The temperature dependence of the ct-tt' polarised in- 
tensities are shown in Fig. [TWa) where they are seen to 
disappear above T;v2=28 K. Magnetic satellites were also 
identified at these wave- vectors in the neutron scattering 
investigation by Kenzelmann et al.^'^ and were attributed 
to an induced moment on the Tb ions, oriented along 
the a axis, from the onset of the Mn magnetic structure 
in the cycloid phase. From Eq. [5] it can be seen that 
the observed azimuthal dependence is consistent with a 
magnetic moment directed along the a axis. Detailed 
analysis described in Sec. [D] substantiates this further, 
but also indicates that the weak pre-edge feature in the 
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unrotated channel shown in Fig. [14] may reflect the con- 
tribution from an order parameter with both odd time 
and inversion parity. 

In summary, in the cycloidal phase additional XRS 
satellites appear at {0,k ± q, 0), k odd, for photon en- 
ergies around the Tb L3 edge, but not the Mn K edge. 
They are found predominantly in the rotated polarization 
channel, and are most probably El — El XRMS from a 
polarisation of the Tb 5d states along the a axis induced 
by the onset of the Mn cycloidal order. These findings 
are in agreement with the model of Kenzelmann et al. 
for this phase. 

It is of course important to reconcile these XRS ob- 
servations at the Tb edge with our failure to observe 
NRXMS at (0, A: ± q, 0) k odd, discussed in Sec. EH 
This may be achieved if we suppose that the Tb 4f mo- 
ments are only weakly polarised in the cycloidal phase. 
This combined with the highger background (Fig. [5^d)) 
observed in the a-a' channel where the peak would be 
expected to occur (Eq. [ij would then act to render the 
peak unobservable. 



6. Low temperature phase: Tb L edge 

We have additionally investigated the XRS below 
Tn3 = 7 K, where the Tb moments order, with the dis- 
tinct wave-vector of (0,qT6,0) with qrb=0.42 b*. How- 
ever, in this phase, energy scans taken at QTb find only a 
weak XRS response, as shown in Fig. [12] (a). The peak 
in the energy scan at the Tb L3 edge is close to the max- 
imum of the fiuorescence spectrum and allows us to asso- 
ciate the XRS as arising from El — El transitions, prob- 
ing the Tb Sd states. For completeness, we also show in 
Fig. [15] (b) an energy scan of q^b taken at the Mn K-edge, 
where no significant XRS can be identified. This implies 
that, within the accuracy of our measurements, we find 
no evidence for an induced polarisation of the Mn ions 
from the onset of the Tb order. 

A scan of the wavevector transfer parallel to [OfcO] 
through q^b provide us with more information about the 
weakness of the XRS, Fig. [TH] (a). This scan reveals that 
the peaks associated with qyb are very broad in recipro- 
cal space, even at the lowest temperature measured of 2 
K and are therefore consistent with only short range or- 
der. Presumably, this effect arises from the competition 
of magnetic order at this position with the induced Tb 
ordering we identify at (0, fc ± q, 0) k odd. Indeed, it is 
intriguing that these XRS investigations appear to find 
that the major polarisation of the Tb valence electrons 
arises from the Mn id order and not from the ordering 
of the Tb 4/ ions. The short ranged order at q^b has 
also been reported in forgoing neutron investigation o^^'^^ . 
It may be of significance that the coupling of Mn and 
Tb moments have been invoked as higher order terms 
in the nearest-neighbour and next-nearest-neighbor(NN- 
NNN) model to explain the emergence of the incommen- 
surate spin structure in the heavy rare-earth ReMnOa 
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FIG. 15: (Color online) (a) Energy scan taken at the Tb 
ordering wave-vector qTi,=0.42 at T=2K in the a — n' channel, 
where only weak XRS is observed, (b) Energy scan taken at 
the Mn ordering wave-vector qAfn=0.28 at T=2K, where no 
significant XRS is observed. 

systemai^. 

The azimuthal dependence of (0, 4 — q^b, 0) is shown in 
Fig. lllf c). and has identical azimuthal symmetry as those 
recorded at (0, fc±q, 0) k odd as shown in Fig.[31[b). This 
comparison indicates that the low temperature ordered 
Tb magnetic moments are orientated along the a axis. 
The temperature dependence of the intensity of (0, 4 — 
qTbjO) is shown in Fig. [TSJ^b), where it disappears above 
T/vs in agreement with neutron scattering results. 



7. Charge satellites 

We have also investigated the satellites which occur at 
2qMra in TbMnOs. We find that these satellites around a 
given A- type magnetic wave- vector such as (0, 4 — q, 1) do 
not resonate at either the Tb or Mn absorptions edges. 
Of significance is the fact that we also find 2qMn satellites 
associated with (0,4 — q, 0). This result suggest that 
the charge degrees of freedom also couple to the lattice 
causing the charge order satellites that we observe and 
indicating an intricate interplay of magnetic and charge 
lattice coupling in this material. We also found that these 
2qAfn peaks do not resonate at either absorption edge. 
Key data on the wave- vector scans are shown in Fig. [T7] 
for the resonant and non-resonant 2qjv/n satellites. 
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FIG. 16: (Color online) (a) k-scan taken at the Tb ordering 
wave-vector qT6=0.42 at T—2K in the a — n' channel. The 
broad width of the scattering indicates that underlying mag- 
netic order is only short ranged ordered at this temperature, 
(b) k-scan taken at the Mn ordering wave-vector qMn~0.28 
at T=2K indicating that the induced polarisation on the Tb 
ions are longed ranged order at this position. 



It is also noteworthy that we do not find any evidence 
for charge ordered satellites associated with the C-type 
satellites such as (0, 3 — q, 0). 



D. Analysis and Discussion 

To obtain a deeper understanding of our resonant scat- 
tering investigations, wc have performed ab initio elec- 
tronic structure calculations using the FDMNES cod^^ 
which allowed us to evaluate the XRS spectra, including 
possible contributions from E1~E1, E1—E2 and E2—E2 
processes. 

It is important to point out that El — El XRS can, in 
most cases, be easily identified because the energy max- 
imum occurs very close to that in the fluorescence spec- 
trum. However, XRS features which occur at energies 
corresponding to the pre-edge threshold may arise from 
El - E1,E2 - E2 and El - E2 events and cannot be 
simply interpreted from their position in energy. 



1. FDMNES code 

FDMNES is a software package that allows the cal- 
culation of the intensity of the diffracted reflections in the 
vicinity of absorption edges of the elements present in the 
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FIG. 17: (Color online) Summary of scans of the wavevec- 
tor transfer along [0,fc,0]. (a) and (b) XRS data around the 
primary satellite positions (0, fc ± qA/n, 0) at the Tb L3 edge 
for azimuthal angles of tp= and 90° in the ct — tt' channel, 
(c) Non-resonant, charge satellites at (0, fc ± 2qMn,0) in the 
a ~ a' channel. 



material of interest. The code may use various methods 
to calculate the excited states probed during the reso- 
nant process and in our calculations, we employed the 
multiple scattering theory. The final states are calcu- 
lated in a fully relativistic way, including the spin-orbit 
contribution. Then the matrix elements governing the 
transition process between the initial and final states are 
calculated considering the polarization conditions of the 
incoming and outgoing photon. This photon electric field 
is expanded in order to calculate both the electric dipole 
{El) and electric quadrupole {E2) contributions and con- 
sequently for the resonant process the El — El, El — E2 
and E2 — E2 components. At the end the structure fac- 
tors are calculated and an expansion in a spherical ba- 
sis is performed, so at the various diffraction satellites 
the different contributions coming from all the order of 
scattering can be separately provided. In this way one 
obtains the monopole (F'^), dipole (F-^) and quadrupole 
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FIG. 18: (Color online) Azimuthal dependence of XRS at 
the (0,3-qMn,0) satellite in the vicinity of the Tb edge 
recorded on ID20. The a — n data were taken at an energy 
corresponding to an El — El event around 7.520 keV, while 
the a — a data were taken 7 eV below this energy. The lines 
have been calculated using the FDMNES package. In the 
unrotated channel, the main contribution is calculated to be 
from El — El XRMS associated with a splitting of the Tb 5d 
bands induced by the cycloidal order on the Mn sublattice. 
For the unrotated channel the weak pre-edge peak is calcu- 
lated to be El — E2, arising from an anapole, i.e. a multipole 
which is odd with respect to both time and parity. The data 
in the unrotated channel have been multiplied by a factor of 
five. 



(F^) El — El contributions, the time even (+) and time 
odd (— ) dipole (F^±), quadrupole (f and octupole 
(F^±) El — E2 contributions and the monopole, dipole, 
quadrupole, octupole and hexadecapole E2 — E2 contri- 
butions. For the El — E2 component the non magnetic 
and magnetic parts can be given separately, so it is pos- 
sible to precisely deduce which terms in particular con- 
tribute to the XRS intensity. 

In practice one has to provide to the code the crystal- 
lographic structure and a guess of the electronic struc- 
ture, including the spin polarization. From this the Dirac 
equation is solved. At low temperature neutron diffrac- 
tion experiments have shown that the crystallographic 
structure does not change by very much2^ from the high 
temperature Pbnm structure. The main problem in the 
context of the present study is that even minute changes 
in crystal structure greatly effect the XRS. We chose to 
keep the Pbnm structure we built a supercell, four times 
bigger along the b direction to mimic the magnetic mod- 
ulated structure of Kenzelmann et al. In this way the 
corresponding new diffraction vectors (/i, k ± q,l) with q 
= 0.25 b* is not far from the incommensurate one mea- 
sured. In order to model the data in the cycloidal phase 
below we kept the same diffraction vector with q 



= 0.25 b* and allowed the magnetic moment to modu- 
late in the b — c plane as described by Kenzelmann et 
al^. Making these approximations we used FDMNES 
to calculate the XRS intensities at the wave-vectors cor- 
responding to A-type (0,4 ± q, 1), i^-type (0,4 ± q, 0) 
and C-type (0, 3 ± q, 0) order, which were extensively in- 
vestigated in our experiments. In all of the following 
discussion, intensities below a threshold of 1 x lO^^rg 
were taken to be zero. This is a reasonable approxima- 
tion as from the FDMNES calculations, the intensity 
at (0, 4±q, 1) at the Mn K-edge is 5 x IQ-'^r^. The corre- 
sponding intensity measured in our experiments at ID20 
at this wave- vector was about 100 cts/second, so that in- 
tensities of the order 500 times below this are beyond the 
sensitivity of our measurements. For the coUinear phase, 
the calculated intensities of all possible multipolar con- 
tributions were found to be negligibly small at the Tb L3 
edge, while the results for the Mn K edge in this phase 
are given in Tab. [D Corresponding calculated intensities 
for the cycloidal phase are given in Tab. HIl and [1111 The 
El — El intensities are given 10 eV above the Fermi level 
(0 eV), while the El - E2 and E2 - E2 intensities are 
given at a pre-edge value taken 5 eV below the Fermi 
level. 



2. CoUinear phase 

We first consider our results for the coUinear phase, 
Tn2 < T < T/vi. Here our XRS data reveals two sets 
of satellites: ^-type, such as {0,k ± q, 1); and i^-type 
(0, fc±q, 0), both with k even. These two sets of satellites 
showed resonant behaviour at both the Mn K and the 
Tb L3 edges. The results of the FDMNES calculation 
at the Mn if -edge, using the b-axis modulated magnetic 
structure of Kenzelmann et al.^'' and the Pbnm structure, 
are shown in Fig. [TWa)-(c) for the ^-type, F-type and 
C-type structures, respectively, while the decomposition 
of the intensity into its multipolar components is given 
in Tab. HI The calculations predict finite El — El F^ at 
the A- type wavevectors, of the order 5 x lO^'^rg, and that 
the E2 — E2 intensities are negligibly small. The pre-edge 
peaks shown in Fig. [19] (a), are all El — El in origin. The 
El — E2 multipoles are negligible in these calculations, 
since the Mn ions occupy sites with inversion symmetry 
in the Pbnm structure. Interestingly, the code also pre- 
dicts some finite El — El XRMS arising from the F^ 
multipole at the F-type (0, 4 ± q, 0) wave- vector, which 
is an order of magnitude weaker than that calculated at 
(0,4 ± q, 1). This result is in contradiction with our x- 
ray experiments, where we find similar XRS intensities 
for the A- and F-type peaks at the Mn K-edge as shown 
in Fig. [Sj (b) and (d), but adds weight to the idea idea 
that the XRS observed at the F-type peaks is not pre- 
dominantly magnetic in origin. Shifting the Mn ions from 
their nominal positions in the calculations had the effect 
of causing larger scattering at the F-type wave-vector, 
but yielded azimuthal dependences inconsistent with our 
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FIG. 19: (Color online) Representative results of FDMNES calculations for TbMnOs in the vicinity of the Mn K edge for 
the coUinear ((a)-(c)) and cycloidal ((d)-(f)) phases. The calculated x-ray absorption near edge structure (XANES) is given by 
the dot-dashed line in (a). The figure allows for a comparison of the calculated intensities of the A- type ((a) and (d)), -F-type 
((b) and (e)) and C-type ((c) and (f)) studied in our experiments. The cr — tt' and o — a' channels are represented by the solid 
red and dashed blue lines, respectively. 



experimental data. The FDMNES calculations predict 
no significant XRS at the C-type peaks as shown in Fig. 
[12] (c), in agreement with our experimental observations 
in the coUinear phase. 

Calculations using FDMNES of the XRS in the 
coUinear phase at Tb L3 edge were also performed. The 
calculated XRS spectra are shown in Fig.[20](a)-(c) where 
they are seen to be negligibly small in the sense described 
above. Thus the calculations do not reproduce our exper- 
imental results, where we observed significant XRS inten- 
sities at both the A- and i^-type peaks in the coUinear 
phase at the Tb L3 edge. A possible explanation for 
this discrepancy is that the FDMNES calculations do 
not predict a magnetic polarisation on the Tb 5c? states 
induced from the Mn coUinear magnetic structure. How- 
ever, it is noteworthy that the calculations do identify 
a polarisation on the Tb ions in the cycloid phase (dis- 
cussed below) for the C-type reflections, induced by the 
Mn cycloid magnetic structure. This induced polarisa- 
tion is in agreement with our experimental results and 
the neutron scattering results of Kenzelmann et al. In 
the light of this fact, another plausible scenario is that 



the XRS observed at the Tb edge for both the A- and 
-F-type peaks are not predominantly magnetic in origin. 
Instead they might arise from a charge multipole order 
parameter, which do not appear with measurable inten- 
sity in the calculations because of assumptions made in 
the model, such as the restriction to the Pbnm space 
group. Neither experiment nor calculations yet provide 
a conclusive answer to this issue. 



3. Cycloidal phase 

For the cycloidal phase at the Mn K edge, the cal- 
culated XRS spectra are shown in Fig. [TWd)-(f) and 
the corresponding multipole intensities are summarised 
in Tab. [Til These calculation predict significant XRMS 
intensity (1 x lO^'^ r§) for the A-type and weaker inten- 
sities for the F-type (5 x 10"^ Tq) and C-type (1 x 10^^ 
Tq) peaks, all of which are magnetic dipole in origin, 
El — El F^. As with the comparison in the coUinear 
phase, there is considerable disagreement between exper- 
iment and calculation in the ratio of intensities at the 
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A- and F-type peaks. Although the calculations indi- 
cate that the i^-type intensities are around two orders of 
magnitude smaller, our experiments find similar XRS in- 
tensities for the two peaks. Again, this may be an indica- 
tion that the XRS we observe at the F-type peaks is not 
magnetic in origin. The calculations predict a very weak 
XRS at the C-type wave-vectors, in disagreement with 
our experimental results where these peaks were found 
to be absent. However, this weak intensity may be at 
the limit of our experimental sensitivity. Again, as in the 
collinear phase at the Mn K edge, all El — E2 XRS is 
calculated to be negligibly small at this absorption edge 
because the Mn ions are assumed to occupy positions of 
centre inversion symmetry in the Pbmn space group. 

For the ferroelectric phase, T/v3 < T < Tn2, the most 
striking feature of our XRS data is the appearance of 
C-type, (0, ± q, 0), k odd, satellites for photon ener- 
gies around the Tb L3 edge. In this phase, the XRS 
cross-section at the Tb L3 was calculated by inputing 
the cycloidal magnetic structure. Here the FDMNES 
code predicts a measurable scattering intensity due an 
El - El XRMS process in the a-n' channel of 4 x 10""*. 
Remarkably the approximately sin^ tp azimuthal depen- 
dence predicted in this channel appears to account for the 
data shown in Fig.fm^b). The FDMNES code actually 
reveals that scattering at this wavevector in the rotated 
channel arises from a small magnetic net polarization of 
the Tb 5d bands along the a-axis induced by the on- 
set of cycloidal magnetic order on the Mn sublattice. In 
contrast, only weak XRS intensities were calculated at 
corresponding wave-vectors at the Mn K edge. 

The azimuthal symmetry at (0, fc±q, 0), k odd, is actu- 
ally more complex than the simple sin^ form suggested 
by the data in Fig. [TlT b) for the a — tt' channel. A much 
more extensive set of measurements were performed with 
the results shown in Fig. [181 Not only is the azimuthal 
dependence clearly more complex, with one-fold rather 
than two-fold symmetry, but in addition these experi- 
ments revealed the existence of a weak component to the 
scattering in the unrotated {a — a') channel. Most im- 
portantly, this component has an azimuthal dependence 
that peaks in anti-phase with the rotated (ct — tt') compo- 
nent (Fig. [T8)) . and resonates 7 eV below the position of 
the El — El dipole resonance (Fig. [T3]). The latter fact 
indicates that it arises either from an E2 — E2 or from an 
E1 — E2 process. Calculations using the FDM NES code 
indicate that, although the E2 — E2 and El — E2 spectra 
have similar azimuthal dependencies, the a — a' signal ac- 
tually arises mostly from an E1—E2 event at this position 
(Tab. |TT|. Specifically, this component is consistent with 
the existence of a finite F^^ term, a time-odd, parity- 
odd multipole, which is also known as the anapole. The 
calculations indicate that the E2 — E2 contributions are 
very small at this position because the FDM NES band 
structure calculations predict that the induced splitting 
is largly in the Sc? and not the 4/ band. 

At this stage it is not possible to make a definitive 
conclusion regarding the observation of an anapolar con- 



tribution to the XRS in multiferroic TbMnOa. For one 
thing, the agreement between the data and the calcula- 
tions based on the FDMNES code are qualitative only. 
While the code does a good job of accounting for the az- 
imuthal dependences (Fig. [T8)) and spectra (Fig. [14] and 
Fig. [20jg)-(h)), it overestimates the relative intensity of 
the (ct — cr') by a factor of roughly 5. However, it is prob- 
ably reasonable to assert that, given the complex nature 
of the calculations, that the agreement is about as good 
as might be achievable. 

We reiterate here that the FDMNES calculations pre- 
dicting an anapolar contribution at the C-type reflections 
assume a Pbnm space group. For this structure the Tb 
ions lack a centre of inversion and even in the absence 
of any ferroelectric order an anapole moment on the Tb 
sites is allowed by symmetry. In the FDMNES calcula- 
tions the anapole contribution presumably is induced by 
the non-coUinearity of the Mn moments in the cycloidal 
phase. Whether additional mechanisims contribute to 
the possible anapole moment, such as small shifts of the 
Tb positions, remains an open question. 

Another explanation for the discrepancy between the 
calculated and experimental polarisation dependence at 
the C-type reflections could be that, as we have already 
evoked for the A and F reflections observed at the Tb 
edge in the collinear phase: the C-type peaks might also 
have a charge multipole contribution in the cycloid phase. 
This could have the effect of increasing our measured 
El — El XRS and thus making the experimental polari- 
sation ratio larger than the calculated one. In this respect 
we are unable to draw a definitive conclusion. We also 
note a further weakness with the analysis. At present the 
FDMNES code can only deal with one magnetic ion at 
a time. In addition to the cycloidal magnetic order on 
the Mn sublattice, the model of Kenzelmann et al^ has 
1.0 fiB on the Tb sublattice polarized along the a-axis. 
Including just this latter component in the FDMNES 
code revealed unsurprisingly that at (0, fc ± q, 0), k odd, 
the E2 — E2 contribution is larger than the one from 
El — E2. However, this also introduced significant in- 
tensity at (0, 4±q, 0) with an incorrect azimuthal depen- 
dence and we therefore consider it to be less consistent 
with our experimental data. 



E. Summary 

We have carried out a comprehensive series of x-ray 
scattering experiments on multiferroic TbMnOs using 
both non-resonant and resonant techniques. In conceiv- 
ing of the experiments our hope was to be able to ex- 
ploit the exquisite (and unique) ability of XRS to reveal 
and differentiate between possible multipolar order pa- 
rameters. In particular, the transition from the collinear 
to the cyloidal/ferroelectric state in TbMnOa, with the 
concomitant loss of inversion symmetry, should, at some 
level of sensitivity, be refiected by new features appear- 
ing below the transition arising from mixed El — E2 pro- 
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FIG. 20: (Color online) Representative results of FDMNES calculations for TbMnOa in the vicinity of the Tb L3 edge for 
the coUinear ((a)-(d)) and cycloidal ((e)-(h)) phases. The calculated x-ray absorption near edge structure (XANES) is given by 
the dot-dashed line in (a). The figure allows for a comparison of the calculated intensities of the A- type ((a) and (e)), f-type 
((b) and (f)) and C-type ((c), (d), (g) and (h)) studied in our experiments. The a — -k' and o ~ a' channels are represented by 
the solid red and dashed blue lines, respectively. 



cesses. Evidence for such changes have indeed been found 
at both the Mn K (Fig. [TS]) and Tb L (Fig. ^ edges. 
Hovifever, the job of relating these observed changes to a 
microscopic description of the active multipoles is a con- 
siderable challenge, and one that we have only gone part 
way to completing here. Our hope is that our data will 
serve to stimulate further theory and experiments in this 
class of materials. 

More specifically, aided by NRXMS, which allowed us 
to characterise the magnetic domain structure within the 
sample volume probed by the x-ray beam, our experi- 
ments appear to have uncovered two previously unknown 
order parameters in this material. The NRXMS inves- 
tigation has determined unambigously, that our sample 
has a single A-type magnetic domain only. Netherthe- 
less, XRS scattering is identified at F-type and C-type 
wavevectors, neither of which can be reconciled with 
XRMS from the nominal Mn magnetic order. Since the 
NRXMS results appear to exclude the El - El i^-type 
XRS observed at the Mn K-edge from being magnetic 
in origin, we propose that it most likely arises from an 
induced charge multipole order parameter. If this is cor- 
rect, then of crucial importance are the charge order 



satellites observed at 2qMn about these F-type satellites, 
which suggest an intricate magnetic and charge order 
coupling in the material. We have attempted to further 
analyse the XRS results by implementing the FDMNES 
code, with varying degrees of success. Perhaps the most 
noteworthy outcome of these calculations is the identifi- 
cation of significant E1 — E2 events in the Tb L3 pre-edge 
of the XRS at C-type wavevectors in the magnetic cy- 
cloidal phase. Although our results at this stage should 
be considered as indicative only, our experiments may 
provide evidence for the existence of a novel type of anap- 
olar oder parameter in the rare-earth manganite class of 
multiferroic compounds. 
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APPENDIX A: FDMNES MULTIPOLE 
CALCULATIONS 
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TABLE I: The results of FDMNES calculations of XRS multipole intensities at the Mn K edge, in units of electrons'^, in the 
coUinear phase. All calculations were undertaken for azimuth tp = Q. The El — E2 multipoles are negligible in the calculations 
because the Mn ions are located at cites with centre inversion symmetry in the Pbnm space group. The even order El — El 
and E2 — E2 and the time-even El — E2 terms are explicitly zero because they cannot arise from our magnetic structure. In 
the table zero means the calculated intensities were less than 1 x 10^® i.e. very small, not zero 
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TABLE IL The results of FDMNES calculations of XRS multipole intensities at the Mn K edge, in units of electrons^, in the 
cycloidal phase. All calculations were undertaken for azimuth ■0 = 0. The El — E2 multipoles are negligible in the calculations 
because the Mn ions are located at cites with centre inversion symmetry in the Pbnm space group. The even order El — El 
and E2 — E2 and the time-even El — E2 terms are explicitly zero because they cannot arise from our magnetic structure. In 
the table zero means the calculated intensities were less than 1 x 10"® i.e. very small, not zero 
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TABLE III: The results of FDMNES calculations of all XRS multipole intensities at the Tb L3 edge, in units of electrons^, in 
the cycloidal phase. The ± for the Fl — F2 multipoles refer to time even and time odd multipoles, respectively. All calculations 
were undertaken for azimuth ^ — 0, except the (0,3 ± q, 0) an' which was calculated at azimuth 1/) = 90°. The even order 
FlFl and E2E2 and the time-even Fl — F2 terms are explicitly zero because they cannot arise from our magnetic structure. 
In the table zero means the calculated intensities were less than 1 x 10"® i.e. very small, not zero 
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